High-order cumulants in the counting statistics of asymmetric quantum dots 
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Measurements of single electron tunneling through a quantum dot using a quantum point contact 
as charge detector have been performed for very long time traces with very large event counts. 
This large statistical basis is used for a detailed examination of the counting statistics for varying 
symmetry of the quantum dot system. From the measured statistics we extract high order cumulants 
describing the distribution. Oscillations of the high order cumulants are observed when varying the 
symmetry. We compare this behavior to the observed oscillation in time dependence and show that 
the variation of both system variables lead to the same kind of oscillating response. 

PACS numbers: 72.70.+m, 73.23.Hk, 73.63.Kv 



Current fluctuations in mesoscopic systems allow to 
obtain information on transport that is not accessible 
from the average current alone^. Measurements of the 
shot noise as a first step beyond the mean have been de- 
ployed successfully to examine correlations in the electron 
transport through a metallic islancP and through semi- 
conducting quantum dot systems^, but the extraction 
of higher moments directly from the current fluctuations 
is a demanding taslP and has not yet been achieved for 
transport through quantum dots. In contrast higher mo- 
ments are naturally accessible in the context of full count- 
ing statistics (FCS)P^. With the use of a quantum point 
contact (QPC) as a non-invasive charge detector with suf- 
ficiently fast time response, FCS became experimentally 
feasible in quantum dot physics^H^. In a measurement 
of the 4th and 5th cumulant complex behavior as a func- 
tion of asymmetry and integration time with local min- 
ima was found 11 . A significantly improved experimen- 
tal technique has made the extraction of very high-order 
cumulants possible, revealing an oscillating behavior as 
function of integration time that strongly increases with 
the order of the moment 1 2 * 13 l Such oscillations are also 
predicted theoretically in quantum optics^ and particle 
physics 14 . The theoretical treatment has shown that these 
oscillations are a universally expected phenomenon for 
most physical systems and that in quantum dots they 
should sho w up prominently as function of the barrier 
asymmetry * 12 * 16 *. 

In this letter we present the experimental verification 
of the predicted oscillations in the higher cumulants as 
function of barrier asymmetry in a quantum dot system. 
To this aim we employ real time single electron count- 
ing with a high bandwidth detector. This enables us 
to measure the full counting statistics for electron trans- 
port through a quantum dot on the /is-timescale, allowing 
us to gather a unprecedented large amount of counting 
events. With this large statistical basis we can exam- 
ine the functional dependence of high order cumulants, 
which describe the statistical properties of the system, on 
the tunneling barrier asymmetry and also on the count- 
ing time. The observed oscillations as function of both 
parameters show large similarities that strongly support 
the universal nature of these oscillations. 
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FIG. 1. (a) Average current through the quantum dot ver- 
sus gate voltage at the in-plane gates Gl and G2. Current is 
measured counting electron by electron. A bias of 1.3 mV is 
applied between source and drain contact of the quantum dot. 
Inset: AFM image of the structure, (b) Measured tunneling 
rates Ys (source), Yd (drain) and the electron transport rate 
with respect to the asymmetry parameter a. The asym- 
metry parameter is modified by changing the quantum dot 
gate configuration along the line shown in (a). 



Our device is based on a GaAs/AlGaAs heterostruc- 
ture containing a two-dimensional electron system 
(2DES) 34 nm below the surface. The electron density is 
p = 4.6- 10 15 m -2 , and the mobility is /i = 64 m 2 /Vs. We 
have used an atomic force microscope (AFM) to define 
the quantum dot (QD) and the quantum point contact 
(QPC) structure by local anodic oxidation (LAO) of the 
surface^EHl the 2DES below the oxidized surface is de- 
pleted and so insulating lines are defined. 

An AFM image of our device is presented in Fig. [I] 
The dark lines depict the insulating oxide barriers writ- 
ten by the AFM. The QPC (left area) is separated from 
the QD structure (right area) by an insulating line. The 
QPC can be electrically tuned using the in-plane gate 
G3. The current through the QPC is amplified by a 
current amplifier and detected in a time-resolved man- 
ner with a sampling rate of 500 kHz. To achieve a high 
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bandwidth charge detection a custom low capacitance 
wiring is used to reduce technical noise caused by the 
input voltage noise of the amplifier. The bandwidth of 
the charge detection is only limited by the current ampli- 
fier bandwidth and slightly exceeds 100 kHz. The QPC 
bias is chose n sufficiently small to avoid back-action on 
the QD 19 20 . The QD is coupled to source and drain 
electrodes via two tunneling barriers which can be sepa- 
rately controlled with gate voltages Vqi and Vq2- These 
gates are also used to tune the number of electrons on the 
QD. The sample is placed in a He 3 refrigerator reaching 
temperatures down to 400 mK. 

Setting the device to a situation where electrons can 
enter and leave the quantum dot results in a fluctuating 
charge on the quantum dot. The electron number on the 
dot is sequentially changing from N to N+l and from 
N+l to N electrons. This leads to a corresponding change 
of the potential at the QPC which has a working point 
on the edge of the first conductance step. Thus the QPC 
acts as a sensitive charge detector as each change in the 
potential due to a changing electron number on the dot 
leads to a current change at the QPC. 

The average current versus Vqi and Vq2 as measured 
by single electron counting is shown in Fig. [T^i. A bias of 
1.3 mV is applied to the quantum dot. The white areas 
without electron tunneling events correspond to Coulomb 
blockade with fixed electron number on the quantum dot. 
In between wide stripes with a significant number of tun- 
neling events are observed due to the applied QD bias 
voltage. The charging energy of the quantum dot is 1.6 
meV, the single particle levelspacing is 0.27 meV. The 
nonzero electron temperature leads to an ambiguous tun- 
neling direction at the edges of these stripes where the 
Fermi- level of a lead is in resonance with the dot level. 

To analyze only unidirectional transport, the system 
is tuned to a situation where the transport state of the 
dot is sufficiently far from the source and drain Fermi 
energies, thus allowing no back-tunneling from the dot 
to the source contact or from drain to the dot (compare 
inset in Fig. |T]o) . The symmetry of the tunneling rates 
are studied using both quantum dot gates to keep the 
transport in the unidirectional regime and concurrently 
to change the asymmetry of the tunneling barriers. The 
line drawn in Fig. [TJi indicates the gate voltages used in 
the following. The resulting tunneling rates are displayed 
in Fig. [TJd with respect to the asymmetry factor a = 
(Ts — Td)/(Ts + Td) with Ts and Tjj the tunneling 
rates between source resp. drain and the quantum dot. 
For smaller Vqi and more positive Vqi (upper left corner 
of Fig. |TJl) the system is rather asymmetric with T$ ~ 
20 kHz and T D w 0.2 kHz. Following the line T s drops 
down to 2.3 kHz and increases to 1.2 kHz. Thus the 
asymmetry factor a is varied from 0.98 to 0.3. The total 
rate of electron transfer = Ts • Td/(Ts -\-Td) reaches 
its highest value = 0.9 kHz at a = 0.3 and drops with 
increasing asymmetry a. 

The detection of individual tunneling events allows us 
to extract the full counting statistics of electron transport 




FIG. 2. Cumulants of order m — 2 ... 15 as function of 
asymmetry factor a. While the low order cumulants show a 
monotonic behavior the cumulants with order m > 4 reveal 
distinct oscillations with factorial increase in amplitude (up to 
10,000 for the 15th cumulant). The cumulants are determined 
for a constant dimensionless counting time interval Ts • to = 
3.2. 

through the quantum dot. Before going into this analysis 
it is useful to introduce a dimensionless time Ts-t which is 
equivalent to setting Ts equal to unity. This is necessary 
to separate time and asymmetry dependence as accord- 
ing to Ref.^time dependent oscillations are universal on 
the timescale Ts-t. Then the system is solely governed by 
the asymmetry parameter a. Thus the statistical prop- 
erties are also completely determined by this single sys- 
tem parameter a. Of course we have to note that in 
our experiment another rate exists, namely the detector 
bandwidth T^et- However, if we fix this rate relatively 
to we again have a system whose statistical proper- 
ties are solely varied by the asymmetry. We realize this 
fixation of the detector bandwidth by a post-selection of 
detected events with an artificial detector bandwidth of 

r d et = 4-r 5 . 

We will now examine the counting statistics of the 
quantum dot. This can be done either by direct anal- 
ysis of the whole probability distribution of the number 
of transferred electrons n in a certain time interval or 
more common by the cumulants c m of this distribution. 
The first cumulant is the mean of the number of trans- 
ferred electrons, c\ = (n), the second, C2, is the variance, 
and the third is the skewness. Cumulants of higher or- 
der m are very sensitive to the details of the transport 
process. Recently it was found experimentally that the 
high-order cumulants of quantum dot transport oscillate 
as function of the counting time interval!^ In the same 
paper it is stated that such oscillations are expected uni- 
versally for nearly any system when changing a relevant 
parameter, e.g. the symmetry of the tunneling rates in 
a quantum dot system^. For our experimental study of 
this predicted symmetry dependence of high-order cumu- 
lants we varied the asymmetry parameter a and measured 
for each a about 1 million tunneling events. This statis- 
tical sample enables us to determine cumulants up to the 
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FIG. 3. Time and symmetry dependence for the 15th cumu- 
lant. Oscillations can be observed for varying time base Vs-to 
and symmetry a. Two profile lines are shown displaying the 
substantial similarity of both types of oscillations. 

15th order with high precision. 

In Fig. [2] the normalized cumulants d^/c 1 are shown 
as function of the asymmetry factor a for m = 2 ... 15. 
For high asymmetry (a — » 1) all nomalized cumulants 
start from a value of 1 due to the Poissonian limit of a 
single dominant barrier. The first cumulants then drop 
monotonically towards the symmetric limit. For m > 4 
a more complicated behavior can be observed. Begin- 
ning with the fourth order cumulant a minimum starts 
to develop, shifting to higher asymmetry and increasing 
in amplitude with rising order. For the 6th order an addi- 
tional maximum appears that also gets more prominent 
with rising order of the cumulants. This trend continues 
for even higher cumulants, the number of oscillation in- 
creases and also the amplitude rises dramatically. For the 
12th order values up to 145 are reached and the ampli- 
tude of the 15th cumulant already exceeds 10000. This 
factorial increase was shown already for oscillations as 
function of time^ and appears here again, showing the 
universal nature of the oscillations. 

The asymmetry dependence shown in Fig. [2] was evalu- 



ated for a certain choice of the dimensionless length of the 
counting interval. We can also determine the dependence 
on the dimensionless counting interval for each value of 
the asymmetry a, thereby mapping out the cumulants as 
function of both parameters. The result obtained for the 
15th cumulant C15 is shown in Fig. [3^i. The figure nicely 
reveals the similarity in the dependencies on asymme- 
try and on dimensionless time - the general features are 
nearly symmetric with respect to the diagonal. For a 
further illustration of these similarities we plot two pro- 
files in Fig. [3}d and taken at fixed asymmetry resp. 
fixed normalized time. Both curves, C15 as function of 
time Vs-to resp. asymmetry a, show oscillations of very 
similar shape and amplitude. This correspondence be- 
tween the different parameters nicely demonstrates that 
the occurrence and the strength of these oscillations are 
of universal nature^. It is also interesting to examine the 
behavior at long times or towards the case of symmetric 
rates, a = 0. For long times the dependence on Ts • to 
becomes weak and in the long time limit the cumulants 
will become sole function of the asymmetry. Towards 
a = the dependence on asymmetry flattens out and 
the cumulants depend mainly on time. This transitions 
again reveal the similarities in the dependence on these 
different parameters and support the underlying univer- 
sal nature of the oscillatory behavior. 

In conclusion, we have realized a high-bandwidth de- 
tection of single electron transport in a quantum dot with 
a sufficiently large event number to examine the symme- 
try dependence of the cumulants of the counting proba- 
bility distribution up to the 15th order. We have verified 
the theoretically predicted oscillation of the high order 
cumulants as function of the quantum dot barrier asym- 
metry. Furthermore we have performed a 2d-parameter 
study of the cumulants as function of both asymmetry 
and time. Our experimental results nicely reveal the sim- 
ilarity of the dependence on these different parameters, 
thereby demonstrating the underlying universal nature 
of the cumulant oscillations. 
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